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ABSTRACT 

We present the first dynamical mass estimates and scaling relations for a sample of Sunyaev-Zel'dovich effect 
(SZE) selected galaxy clusters. The sample consists of 16 massive clusters detected with the Atacama Cosmol- 
ogy Telescope (ACT) over a 455 deg 2 area of the southern sky. Deep multi-object spectroscopic observations 
were taken to secure intermediate-resolution (R ~ 700-800) spectra and redshifts for w 60 member galaxies on 
average per cluster. The dynamical masses Mitxic of the clusters have been calculated using simulation-based 
scaling relations between velocity dispersion and mass. The sample has a median redshift z = 0.50 and a median 
mass M-2/ooc — 1 1 X lO 14 /i7 ( )M0 with a lower limit M200C — 5 x 1O 14 /i7oM , consistent with the expectations for 
the ACT southern sky survey. These masses are compared to the ACT SZE properties of the sample, specif- 
ically, the central SZE amplitude yo, the Compton signal within a 0.5' pixel yo.'5> an d the integrated Compton 
signal I2000 which we use to derive SZE-Mass scaling relations. All SZE estimators correlate with dynamical 
mass with low intrinsic scatter (11% - 16%), in agreement with numerical simulations. The influence of dy- 
namically disturbed clusters on these scaling relations is also considered. Using the 3-dimensional information 
available, we divide the sample into relaxed and disturbed clusters and find that ~ 50% of the clusters are dis- 
turbed. We conclude that disturbed systems do not significantly bias the scaling relations but might modestly 
boost their scatter. 

Subject headings: galaxy clusters: general, redshifts — galaxy clusters: dynamics, masses — SZ effect: scaling 
relations 



1. introduction 

t Based in part on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, under pro- 
grams 084. A-0577 and 086. A-0425. 

it Based in part on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy, Inc., under a cooperative agreement with the NSF on behalf of 
the Gemini partnership: the National Science Foundation (United States), 
the Science and Technology Facilities Council (United Kingdom), the Na- 
tional Research Council (Canada), CONICYT (Chile), the Australian Re- 
search Council (Australia), Ministerio da Ciencia e Tecnologia (Brazil) and 
Ministerio de Ciencia, Tecnologia e Innovation Productiva (Argentina). 

t Visiting astronomer, Gemini South Observatory. 



Studies of clusters of galaxies have had a wide impact on 
our unders tanding of galaxy formation and cosmology (see 
|Voit||2005| for a review). They are a unique laboratory for 
studying the effects of the environment (high densi ty, gas 
pressure, coll i sions, etc.) on galaxy evolution (|Butcher &| 
|Oemler|[T984l |Balogh et al. 1999; [Hansen et al.||2009K ~At 
the same time, number counts of galaxy clusters, sensitive to 
the amplitude of matter fluctuations, can provide constraints 
on various cosmological parameters (|Bahcall & Fan]|1998 



Evrard et al 12002} |Vikhlinin et al.|2009HMantz et al.|20i0b|c 
Rozo et al.J|2010). An accurate determination of the latter 



requires that we know the mass and redshift distributions of 
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clusters with good precision. 



Th e Sunyaev-Zel'dovich Effect (SZE; Zel'dovich & Sun- 
Sunyaev & Zel'dovich 1970|l is a distortion in 



yaev 



1969 



the Cosmic Microwave Background (CMB) temperature pro 
duced by inverse-Compton scattering of CMB photons as 
they interact with the hot electrons of the intracluster medium 
(ICM) of a galaxy cluster. Its surface brightness is indepen- 
dent of redshift, and its strength is proportional to the line-of- 
sight column density times the electron temperature. The SZE 
is a powerful tool for detecting massive clusters to high red- 
shifts (see, e.g., the reviews by Birkinshaw 1999; |Carlstrom| 
|et al.|2002) . 

Early measurements of the SZE were achieved with tar- 
geted observations of known clusters. These revealed the 
power of SZE studi es, reaching from gas physics and in ner 
structure of clusters (Greg oet al.|2001||Benson et al.|20 04), to 
cosmological parameters such as the Hubble constant (Birkin- 
|shaw et al.|1991|[Hughes & Birkinshaw 199 8]) and the energy 



density of matter in the universe, Um (Grego et al. 2001 ). 
Large SZE surveys over cosmologically significant areas of 
the sky have recently com e to fruition as the Atacama Cos- 
mology Telescope (ACT, |Fowler et al.||200"7} |Swetz et al. 



201 1\ and the South Pole Telescope (SPT, ICarlstrom et al. 



2009 ) have begun scanning large areas of the sky at millime- 
ter wavelengths. The Planck satellite (Taub er et al.|20 10) is 
conducting an all-sky survey and has recently released the first 
all-sky sample of SZE-selected galaxy clusters (Planck Col- 



|laboration|20 1 1 a[ ). T he first cluster detec tion s with ACT and 
SP T are p resented in |Hincks et al.| ( |2010| ) and |Staniszewski et] 
[aE]d2009), respectively! 

The rapidly growing SZE cluster samples have the poten- 
tial t o place strong constrain ts on cosmological parameters 
(e.g., Battye & Weller 12003). Both numerical simulations 



qSpringel et al.|2001a||da Silva et al.|2"004l |Motl et al.|2005[ 
Nagai|2006[|Battaglia et al.|201 1) and analytical stu dies ( |Reid| 
& Spergel|2006t|Ashfordi|2008[ |Shaw et al.|2008) suggest a 
tight correlation between cluster mass and SZE signal. On 
the other hand, biased mass estimates can have a la rge im- 
pact on cos mological parameter determination (e.g., |Francis| 
|et al.|2005). By limiting thei r study to the high-significance 
clusters, |Sehgal et al.| ( |201 1[ ) have shown the power of the 
ACT sample in constraining cosmological parameters, par- 
ticularly the dark energy equation-of-state parameter w and 
the root-mean-square (R MS) mass fluctuations on a scale of 
8/r'Mpc, ag. Likewise, Vanderlinde et al. I (120101 have used 
SPT data to set cosmological constraints, with similar find- 
ings. They have also shown that these improvements can be 
achieved only in the presence of a well-calibrated scaling re- 
lation between mass and SZE signal. To assess the scaling of 
SZE signal with mass, independent means of measuring the 

mass are crucial . 

Ben son et al.| ( |201 1) used X-ray observations in combina- 
tion with SZE measurements to derive an empirical scaling 
relation between mass and SZE signal. This allowed them 
to confirm that SZE-selected samples of clusters yield sig- 
nificant improvements when added to other datasets to con- 
strain cosmological parameters. While X-ray observations 
have proven to be an effective way of measuring cluster 
masses , and have been exploited to characterize the SZ E sig- 
nal (|LaRoque et al.|2006[|Bonamente et al.|2008||Andersson| 
|et al.||2011[|Merin et al.|2011||Planck Coflaboration||2011b) , 
they are do not provide truly independent mass estimates from 
SZE measurements, since both rely on the properties of the 



gas in the ICM and should be affected by similar physics. 

The velocity dispersion of cluster member galaxies is one of 
the most widely used methods for constraining cluster mass, 
and is independent of the properties of the gas in the ICM. It 
takes into account the galaxy distribution and relies, to some 
extent, on the assumption that the clusters are relaxed (i.e., 
virialized). Until recently, however, mass measurements to 
independently calibrate t he SZE signal with mass have come 
from optical richness (Menanteau & Hughes 2009; High et al. 
|20T0l |Menanteau et al.||2010a| |Pla nck Collaboration 2011c) 
and len sing analyses ( Sea lfon et al.|2006[|Umetsu et al.|2009 
Marrone et al.||2011|). | Hand et al.| ( |2011| ) presented stackec 



ACT data in the directions of Luminous Red Galaxies (LRGs) 
from t he Sloan Digital Sky Survey (SDSS) Data Release 7 
(DR7, Abazajia n et al.|2009 1 using optical luminosity-based 
masses. This approach allowed them to probe the SZE signal 
fr om lower mass systems than otherwise possible. 

|Rines et aL] ([2010 ) presented the first statistical comparison 
between dynamically-estimated masses and integrated SZE 
signal from a sample of 15 nearby (z < 0.3) galaxy clusters, 
showing that masses thus determined and the integrated SZE 
flux are correlated at the s» 99% confidence level. Further- 
more, they estimate that the significance is higher than that of 
the correlation between SZE and weak lensing masses from 



Marrone et al. ( 2009 ), probably because of the small apertures 
used in the latter study. However, since their sample was not 



homogeneously selected, |Rines etaL ( 2010| ) do not account 
for observational biases in their sample and do not report a 
formal scaling relation between mass and SZE flux. 

In this work we present spectroscopic redshifts and — for 
the first time — dynamical masses of a sample of clusters of 
galaxies selected with the SZE. These clusters were o bserved 
by ACT in i ts 2008 southern sky survey at 148 GHz (Marriage 
et al.|20l"Ta} , and were optically confirmed by [Menanteau et 



al.| ( |2010bj ). We use a variety of SZE diagnostics to assess the 
scaling with dynamical mass and thus present the first robust 
scaling relations between dynamical masses and SZE signal 
for a sample of SZE selected clusters. 

Throughout this work we us e a flat ACDM cosmology 
consistent with WMAP-7 data ( |Komatsu et al.pqil) , with 
D, A = 0.73, tt M = 0.27 and H = 70h 70 km s" 1 Mpc" 1 . Masses 
and integrated SZE signals are estimated within a radius r200c 
which encloses a density 200 times the critical density of the 
Universe at the redshift of the cluster, p c {z) = 3i/ 2 (z)/87rG. 
All quoted errors are 68% confidence intervals unless other- 
wise stated. 

2. OBSERVATIONS 

2.1. ACT SZE Observations 

ACT is a six-meter off-axis Gregorian telescope operating 
at an altitude of 5200 m in the Atacama Desert in Chile, de- 
signed to observe the CMB at arcminute-scale resolution. It 
has three 1024-element arrays of transition edge sensors oper- 
ating at 148, 218 and 277 GHz. ACT began observing in late 
2007, and su rveyed two regions in the sky with a total area 
of 755 d eg 2 ( |Marriage et atp011b| [Hasselfield etaL|2012, 
|in prep.| ). The processes o f cluster detection and e xtraction 
are thoroughly described in Marri age et al.| ( |201 la| and ref- 
erences therein). In short, the maps are match-filtered and 
convolved with a beta-model profile with j3 = 0.86 with vary- 
ing core radius 9 C from Q'.25 to 4.'0. Cluster signal-to-noise 
ratio (S/N) is measured as the maximum S/N from this set of 
filtered maps. 



ACT Dynamical Masses of Clusters 
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We report on a large spectroscopic follow-up campaign of 
an ACT 148 GHz cluster sample, which was obtained from a 
455 sq. deg. survey of the southern sky. The survey is roughly 
bounded by right ascensions 00 h 12 m and 07 h 08 m and decli- 
nations -56° 1 1' and -49°00'. For further details on the ACT 
observations, data reduction, and cluster detection procedure, 
|Fowler etlE1 ( [2010l ); |Marriage et al.| ( |201 la) . 



Table 1 

ACT-SZE Measurements of Clusters 



see 



2.1.1. The Cluster Sample 

In this study we consider a total of 19 clusters, spanning 
a wide range in mass and redshift. We focus, in particular, 
on the subsample of 16 that were detected by ACT through 
their SZ signal. This subsample contains 15 systems th at were 
detected by ACT in the 2008 single-season maps (Marriage 
|et al.|2011a) and confirm ed optically on 4m-class telescopes 
(Menanteau et al. 2010b), plus one additional cluster (ACT- 
CL J0 52l^ 5104) detected in the new ana lysis of multi-season 
maps (Hasselfield et al. 2012, in prep.) . This latter cluster 
was initially targeted for spectroscopic fol low-up based on its 
optical richness alone (Menanteau et al. 2010a). These 16 
clusters were selected based on a redshift cut or z p hot > 0.35 
and were all discovered with the SZE, with the exce ptions 
of ACT-CL J0330-5227 (X- rays, |Werner et al.||2007l ) and 
ACT-CL J052 1-5 1 04 (optical, Menante auet al.|20T0a) . ACT- 
CL J0330 -5227 is located 12' North-East (NE) of Abell 3128 
(z = 0.06; [C olless & Hewett|[T987l |Katgert et al.][T996T ), but 
|Werner et al.| ( |2007) found it to be an unrelated, background 
cluster at z = 0.44 based on the observed energy of the Fe 
K X-ray emission line using XMM-Newton observations and 
the optical spectrum of the Brightest Cluster Galaxy (BCG). 



pOTOl have shown that the observed SZE sig- 
d to the background cluster. Four clusters 



Hinc ks et al 
nal is clearly relate 

were initially report ed b y SP T ( Staniszewski et al. 200 9") |Van- 
derlinde et al. 2010| see j ffi.l) and studied optically by Menan- 



teau et al.| ( |26l0a) . ACT-CL J0546-5345 is the on ly cluster 



with a dynamic al mas s estimate prior to this study ( Brodwin 
et al.|2010| see ^6X9) 



Of the 16 ACT clusters, 10 are newly discovered. Menan- 
|teau et aL] ( |2010b) confirmed them as clusters with a BCG and 
an accompanying red sequence of galaxies and studied their 
X-ray properties from archival ROSAT data for the 15 clusters, 
plus Chandra and/or XMM-Newton data when available. The 
clusters cover the range ~ 1 -30 x 10 44 ergs _1 in X-ray lumi- 
nosity as measured in the 0.1-2.4 keV band. Photometric 



redshifts were estimated for these 15 clusters by Menanteau 
et al. ( 2010b) . The spectroscopic redshift range covered by 
the sample is 0.28 < z < 1.07 with a median redshift z = 0.50. 

Additionally, we report on three optically selected, high- 
richness ga laxy clusters from th e Southern Cosmology Survey 
(SCS; [Menanteau et al. 2010a). These clusters were part of 
our 2009B follow-up observations before the ACT maps were 
available for cluster detection, and were not detected by ACT. 
They are briefly discussed in §6.2| 



2. 1 .2. Cluster SZE Measurements 

To characterize the SZE produced by each cluster (in the 
148 GHz band) we study three different estimators, corre- 
sponding to different map filters and different apertures for 
extracting the SZE signal. The values for the sample are listed 
in Table □ " 

The first estimator was introduced by Sehgal et al. ( 2011) 
and corresponds to the SZE signal of the brightest 0.5 pixel 
in a map match-filtered with a Gaussian profile of fixed 2' 



Cluster 




z 


)'0'5 3 


>'o 


l200r 








x 10 -4 


x 10~ 4 


xlO" 10 


ACT-CL JO 102- 


-4915 


0.870 


1.80 ±0.22 


3.50 ±0.43 


6.38 ±0.76 


ACT-CL J0215- 


-5212 


0.480 


0.99 ±0.22 


0.79±0.18 


1.62±0.41 


ACT-CL J0232- 


-5257 


0.556 


0.99 ±0.22 


0.61 ±0.17 


1.23 ±0.30 


ACT-CL J0235- 


-5121 


0.278 


1.03 ±0.22 


0.97 ±0.19 


4.22 ±0.89 


ACT-CL J0237- 


-4939 


0.334 


0.73 ±0.22 


0.93 ±0.26 


4.62± 1.33 


ACT-CL J0304- 


-4921 


0.392 


0.88 ±0.22 


1.51 ±0.31 


4.55 ± 1.08 


ACT-CL J0330- 


-5227 


0.442 


1.39 ±0.22 


1.23 ±0.18 


3.92 ±0.58 


A 1 "T TmzL£ 




n s^o 


U.oo zn u.zz 


1 fid -4- O 99 
1 zn u.zz 


9 01-1-0 AQ 
Z.U1 ztz \J.4y 


ACT-CL J0438- 


-5419 


0.421 


1.54 ±0.22 


1.62±0.13 


4.93 ±0.44 


ACT-CL J0509- 


-5341 


0.461 


1.03 ±0.22 


0.41 ±0.13 


0.51 ±0.22 


ACT-CL J0521- 


-5104 


0.675 


0.59 ±0.22 


0.71 ±0.16 


1.23 ±0.29 


ACT-CL J0528- 


-5259 


0.768 


0.73 ±0.22 


0.49 ±0.13 


0.44 ±0.15 


ACT-CL J0546- 


-5345 


1.066 


1.14 ±0.22 


0.91 ±0.14 


1.15 ±0.14 


ACT-CL J0559- 


-5249 


0.609 


0.81 ±0.22 


0.89±0.14 


2.23 ±0.26 


ACT-CL J0616- 


-5227 


0.684 


1.32 ±0.22 


1.00 ±0.15 


2.04 ±0.22 


ACT-CL J0707- 


-5522 


0.296 


0.77 ±0.22 


0.54±0.21 


2.48 ±0.56 



Note. — Redshifts are listed for reference. See Tablef5]for details. 
a Measurements using 2008 ACT data only. Uncert ainties correspond to the 
RMS scatter in ^recovered - }'tnje from simulations; see |Sehgal et al.pHl) . 



FWHM. This estimato r is referred to asywj hereafter. We use 
the values of yof5 from Sehgal ^t al.| ( |201 1) , which have been 
measured using data from the 2008 observations only. 

We also consider the dimensionless, matched filter SZE am- 
plitude, which we denote as yp. The an alysis us ed to deriv e 
these values is presented in Hasselfield et al. (2012, in prep. ); 
it is similar to that performed by |Marriage et al.| (|2011a) 
but now includes the multi-season (2008-2010) ACT observa- 
tions. In brief, the ACT map is filtered optimally to extract the 
amplitude of the temperature decrement of clusters modeled 



by the generalized Navarr o-Frenk- W hite (NFW, Navarro et al. 
|1995) best-fit profile of Arn aud et ak] ( 120 1 0) — hereafter "the 
A10 profile" — with a fixed FWHM of 2'. The correspond- 
ing temperature decrement, ST, is scaled to yo using th e stan- 
dard non-relativi stic SZE frequency dependence (Sun yaev~&| 
Zel'dovich 1980). Our new yo measurements correlate well 
with yo(5, with a Pearson's r-value of 0.75. 

The previous two SZE estimators can be used without any 
prior knowledge of the properties of the clusters to be de- 
tected in a large survey area. Hence, determining the mass 
scaling laws for these estimators is essential for cosmologi- 
cal studies based on the abundance of massive clusters. Once 
the redshift of the cluster is known, an integrated Compton 
signal Y = J ydVl can be measured for each cluster within a 
physical radius. Larger integration areas tend to give mea- 
surements that are more robust ag ainst the effects of cluster 
physics such as AGN f eedback ( |Motl et al.|2005||Nagai|2006l 
Reid & Spe rgel 2006), and to projection effects ( |Shaw et"aT 
2008). Velocity dispersion-mass scalings are usually mea- 



sured at r2ooc (see Eq. fT). The dynamical information there- 
fore yields a measurement of the size of the cluster, which 
is used in the calculation of Y2000 We use the same match- 
filtered map as for the extraction of yo, and the A 10 profile to 
extract each cluster's spherical Y value from within r2ooc- 

2.2. Optical Spectroscopy 
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Table 2 

Spectroscopic Observation Details 



Run 


Semester 


PI 


Tel./Inst. 


Program ID 


Mode 


Grating 


Hours 


Nd 


1 


2009B 


Infante 


VLT7FORS2 


084.A-0577 


Service 


GRIS 3001+11 


14 


3 


2 


2009B 


Barrientos 


Gemini-S/GMOS 


GS-2009B-Q-2 


Service 


R400 G5325 


20 


4 


3 


201 OB 


Infante 


VLT7FORS2 


086.A-0425 


Service 


GRIS 3001+11 


15 


2 


4 


201 OB 


Barrientos/Menanteau a 


Gemini-S/GMOS 


GS-2010B-C-2 


Classical 


R400 G5325 


40 


10 



Note. — N c i is the number of clusters observed in each run. Each cluster was fully observed in one run. 
' Joint Chile/US proposal 



The spectroscopic observations were carried out in two 
semesters, 2009B and 201 0B. Each semester was split into 



two observing runs, one with FORS2 at the VLT ( Appenzeller 



etal.119 98 ) and one with GMOS at Gemini South (Ho oketal. 
2004), both telescopes located in Chile. The details of each 
observing run are listed in Table [2] In total, we had 89 hours 
of observation, during which we collected Multi-Object Spec- 
troscopy for 19 clusters. 

Targets were selected by a two-step process. First, a pho- 
tometric redshift-selected catalog was constructed, includ- 
ing galaxies within ±0.1 of the redshift of the BCG. Within 
this catalog, galaxies were visually selected based on their 
gri colors, with preference given to bright galaxies. All 
our spectroscopic observations cover the wavelength range 
~ 4000-8000A. In this range several spectral features are ob- 
servable at the median photometric redshift of 0.54 (Menan- 



teau et al. 2010b| l. These are mainly the Call K-H absorp- 
tion doublet (at rest-frame wavelength Ao ~ 3950A), which is 
the spectral signature of elliptical galaxies, plus other absorp- 
tion lines such as the G-band (A = 4300A), H/3 (A = 4861 A) 
and the Mgll-triplet (A ( > ~ 5175A), plus the [Oil] emission 
line at rest-frame Ao = 3727 A. The Nal absorption-doublet 
(Ao ~ 5892A) is also observable in the low-z clusters. 

2.2.1. VLT-FORS2 Observations 

The FORS2-2009B observations (Run 1) were aimed at 
newly SZE-detected clusters regarded as "secure" candidates 
detected with A CT in 2008. These clust ers had already been 
reported by Staniszewsk i et al.| ( |2009[ ) as S ZE detections 
and their phys ical properties characterized in Menanteau & 
Hughes (2009). 

Run 3 was mostly focused on getting d etailed informa- 
tion f or ACT-CL J0102-4915 ("El Gordo," penanteau et all 
|201 1) , which was detected as the largest decrement in the 
ACT maps. ACT-CL J0559-5249 was also included in this 
run. 

Runs 1 and 3 were executed in Service Mode in semesters 
2009B and 201 0B, respectively. The instrument setup in both 
runs was the same, using the GRIS 3001+11 grism and I"— 
wide slits, which provides a resolving power R = 660 at A = 
8600A. A total of 18 FORS2/MXU masks were observed for 
the five clusters. Each mask was observed for 40 min, which 
we estimated to be the best compromise between maximizing 
S/N and number of masks. 

FORS2 has a field of view of 6f8 x 6f 8 in the standard reso- 
lution setup, which corresponds to a width of 2517 fcj kpc at 
z = 0.5. 

2.2.2. Gemini-GMOS Observations 

The GMOS-2009B observations (Run 2) were aimed at four 
optically selected clusters from the SCS whose richness-based 



mass estimates sugg ested that they would be d etected by ACT 
in the SZE survey ([Menanteau et al. 2010a[ ). However, as 
mentioned above, only one object was in fact detected by ACT 
(ACT-CL J0521-5104); the other three are discussed in ^2] 
The total integration time per mask was 3600 s (2 x 1800 s). 
Two exposures at slightly different central wavelengths per 
mask were required to cover the two 37-pixel gaps between 
the CCDs which run across the dispersion axis. 

Targets for Run 4 (GMOS-2010B) were selected from the 
sample of clus ters newly discovered by ACT presented in 
Marriage et al. (2011a) and optically confirmed by Menanteau 
|et al.| ( |2010b| ). Run 4 was the only one executed in Classi- 
cal Mode, during five consecutive nights (Dec. 6-10), all with 
clear, photometric conditions and seeing < 0.8". Based on 
our experience in Run 2 we decided to reduce the integration 
time to 2400 s (2 x 1200s) for each mask during Run 4. This, 
coupled with a ~ 20% higher efficiency than Queue Mode, 
allowed us to observe a larger number of masks (and clusters) 
while still obtaining the necessary S/N in the relevant spectral 
lines. 

In both GMOS runs we used the R400_G5325 grating and 
l"-wide slits, providing a resolving power of R ~ 800 with a 
2x2 binning at A — 7000A. In the standard setup GMOS has 
a field of view of 5.'5 x 5'.5 (2036 h^ kpc at z = 0.5). 

2.2.3. Data Reduction 

We have developed reduction pipelines both for the FORS2 
and GMOS data, based on the existing software by ESO and 
Gemini respectively, which work with IRAF/PyRAJff] Cos- 
mic rays are removed using LA. Cosmic ( |van Dokkum|200r| 
with a detection limit of 4.5cr. The wavelength calibrations 
were done using CuAr lamps in the case of GMOS data and 
HeAr lamps for VLT data. The sky is subtracted from each 
spectrum using a constant value determined locally within 
each slitlet. In the case of GMOS data, the individual ex- 
posures are coadded at this point. Finally, the 1 -D spectra are 
extracted from each slit and matched with the input photomet- 
ric catalogs used to generate the masks. 

3. ANALYSIS AND RESULTS 

3.1. Galaxy Redshifts 

Galaxy redshifts are measured by cross-correlating the 
spectra with galaxy spectral templates of the SPSS DR7 us- 



ing the RVSAO/XCSAO package for IRAF dKurtz & Mihk| 
1998); the spectral features in each spectrum have been con- 
firmed with the 2d spectra by visual inspection. We have been 
able to estimate reliable redshifts for ~ 1200 galaxies which 
comprise ~ 80% of all targeted objects. 



The pipeline used to reduce GMOS data — dubbed "pygmos" 
able at http://www.astro.puc.cl/~cjsifon/pygmos/ 



-is avail- 
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Figure 1. The "stacked" result of the shifting gapper method of member 
selection, showing the galaxies in all 16 SZE-detected clusters. The hori- 
zontal axis shows the cluster-centric distance normalized by r2oor for each 
cluster and the vertical axis shows the peculiar velocity of each galaxy, nor- 
malized by the velocity dispersion of the corresponding cluster. Black dots 
show member galaxies, open circles show galaxies rejected by the method, 
and crosses show galaxies with peculiar velocities larger than 4000 kms~'. 
Galaxies with peculiar velocities larger than dSgi are not shown for clarity. 

The median RMS in the wavelength calibration is ~ 0.3A 
and is similar for both instruments. At a central wave- 
length of 6000 A, this corresponds to a velocity uncertainty of 
15kms _1 , which is within the errors of the cross-correlation 
velocity. In particular, the latter is typically A(cz) ~ 40 - 
80km s" 1 , as calculated by RVSAO. It has been established 
experimentally that the true cross-correlation errors are larger 
than those report ed by RVSAO, by a factor ~ 1.7 (e.g., |Quii>| 
tana et aL||2000| l, strengthening the point that the calibration 
errors are well within the velocity measurement errors. 

We have included the me mber catalog for ACT- 
CL J0546-5345 published by |Brodwin "eTaD (|20"l0l> Seven 
galaxies have been observed both by Brod win et al.| l |2010) 
and by us; all redshifts are consistent within 2a. We there- 
fore use our measurements for those galaxies in the following 
analysis. 

3.2. Cluster Redshifts and Velocity Dispersions 

It is of great importance to correctly determine cluster 
membershi p to avoid a biased measurement of the velocity 
dispersion (Be ers et al.||199l) . This is a complicated prob- 
lem and many methods have been developed to handle it. 
In this analysis, membership of galaxies to a cluster is de- 
termined by applying a cut in (rest-frame) velocity space of 
4000 k m s" 1 , and then applying the shifting gapper method 
( Fadd a~et al. 1 1 996) . To do this, we define annular bins around 
the BCG, each of which has at least 15 galaxies and radial 
width > 250/i7 ( l ) kpc. We consider the histogram of velocities 
of member galaxies within each bin. We assume the profile 
is symmetric, and identify the main body of galaxies as those 
wh ose velocity is bounded by gaps of > 500 k ms" 1 . Follow- 
ing Katgert et al.]( |1996| l and Fadda et al. ( 1996), galaxies sep- 
arated from the main body by > 1000 kms -1 are considered 
interlopers and are removed. The selection method is iterated 
until the number of members is stable. This usually happens 
after the second iteration. A total of 948 galaxies (~ 65% of 
all targets) have been selected as cluster members. Most of 



• SOI/SOAR 
o EFOSC/NTT 
o Mosaic/Blanco 



Zspec 

Figure 2. Comparison between spectroscopic redshift s from this work with 
initial gri photometric redshift estimates from Menanteau et al. 12010b}. The 
instruments/telescopes with which each cluster was observed are identified 
in the legend. The dashed line shows z p i,ot = Zspec ■ The dotted horizontal 
line shows the sample selection cut, Zpha = 0.35, and the dotted vertical line 
shows the corresponding z S p ec = 0.35. 

these galaxies show the spectral signatures of elliptical galax- 
ies and do not have emission lines, an d on ly a few emission- 
line galaxies belong to clusters (see §4.3[ ). The galaxies re- 
maining at this point are considered members of the cluster. 
Figure[T]shows the "stacked" result of this method, with mem- 
bers as solid dots. The values have been normalized to allow 
for direct comparison of all clusters. We have explored sys- 
tematic effects coming from the member selection method by 
changing the width of the bins, the number of galaxies per bin, 
and the size of either gap in the shifting gapper. Varying these 
parameters yields results that are consistent with the reported 
velocity dispersions. 

It is clear from Fig. [T] that, while our sampling extends to 
rxyoc in most cases, it is mainly concentrated within the inner 
0.8 r2ooc of the clusters. This is not expected to bias our results 
since velocity dispersion profiles are, in general, flat from 
~ 0.5r20Q C outwards (|Fadda et al.||1996| |Biviano & Girardi] 
2003 1 |Faltenbacher & Diemand|2006|l. We use the biweight 



estimators of location (hereafter z B i) ( [Beers et al.||1990| l for 
the redshift of the cluster and scale, Sbi, for the velocity dis- 
persion. All errors have been estimated with the bootstrap 
resampling technique with 5000 iterations. The redshifts of 
the clusters are presented in Fig. [2) where they are com pared 
to the photometric redshifts of Me nanteau et al.| l f2010b) . The 
median redshift of the sample is z = 0.50. The slightly bi- 
ased photometric redshifts apparent on Fig. [2] are mainly due 
to two factors: the lack of a well-characterized filter response 
function for the telescopes involved in the imaging follow-up 
and the use of only three to four filters for the determination 



o f photometric redshift s ( |Menanteau et al.|2010b I. 

Dane se et al.| ( |1980| l snowed that the observational errors 
on the redshifts of galaxies introduce a bias in the mea- 
sured velocity dispersion. However, for a cluster of ~ 1 x 
10 15 M Q with individual errors as measured in this work (i.e., 
< 100 kms" 1 ), this correction is < 0.1% (and even lower for 
more massive clusters), and it is therefore not considered 
heref] 

5 As mentioned before, the errors calculated by RVSAO are smaller than 
the true cross-correlation errors. Even so, the Danese et al. ( 1 980 1 correction 
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4. DYNAMICAL MASSES 

In this section, we use the velocity dispersions measured 
in the previous section to estimate cluster masses. The dy- 
namical state of each cluster is also studied, including signs 
of substructure and the fraction and influence of emission- 
line galaxies in the cluster population. Both factors can, in 
principle, bias the velocity dispersion and thus the dynami- 
cal mass of the cluster. Moreover, they are not expected to be 
completely independent, since emission-line galaxies are gen- 
erally newly incorporated galaxies, which might mean recent 
(or near-future) mergers involving the main cluster ( Moore et 



|al.|1999[|Book & Benson|20"l0l >. 

4.1. Dynamical Mass Estimates 

The relationship between velocity dispersions and masses 
has bee n the focus of sever al studies. As a first-order ap- 
proach, |Heisler et al.| ( |1985| l studied simple variations of the 
virial theorem and found that they all behave simila rly, and 
that it is no t possible to distinguish among them. Carlberg 
et al. ( T997J compared masses obtained from the virial the- 
orem to those obtained with the Jeans equation in observed 
clusters. They found that the former are biased high by a fac- 
tor of 10-20% and associated this bias with a surface pressure 
correction factor of the same order. 
More recently a nd based on large cosmological simulations 



Evrard et al.| ( |20Q8p concluded that massive (M 2 m c > 10 14 M Q ) 
clusters are, on average, consistent with a virialized state, and 
find a best-fit scaling relation for dark matter halos described 
by NFW profiles in a variety of cosmologies. Accordingly, 
the mass enclosed within rjooc is 



M 2 ooc = 



10 1 



( °DM 



l/c 



M: 



I 



(1) 



where a x5 = 1082.9 ± 4.0 k ms" 1 , a = 0.3361 ±0.0026, 
hioiz) = /j7o\/nA + (l +z) 3 f2 M for a flat cosmology and <tdm 
is the one-dimensional velocity dispersion of the dark matter 
particles within r2ooo which is related to the velocity disper- 
sion of galaxies by a so-called bias factor b v = Sbi/odm- As 
summarized by Evrard et al. (2008), the bias factor as cur- 



rently estimated is (b v ) = 1.00 ±0.05. For consistency with 



previous studies (e.g., Brodwi n" et al.|2010 i, we adopt a value 
b v = 1, meaning that galaxies are unbiased tracers of the mass 
in a cluster. 

The mass values drawn from Eq. [T| are shown in Table [3] 
and the given errors include uncertainties on the cluster recP 
shift, the velocity dispersion, a and CT15. The overall uncer- 
tainty in the mass is dominated by statistical errors which, in 
turn, are dominated by the error in the velocity dispersion. 
The systematics introduced by Eq. [T] contribute < 10% of the 
uncertainties listed in Table [3] The mass from Eq. |T| yields a 
lower value than the virial mass estimator, as Carlberg et al. 
( 1997 1 also anticip ated. 

As indicated by |Evrard et al. ( 2008 ), Eq. [T| holds for pri- 
mary halos, i.e., clusters where a main system ' can be e asily 
identified and substructure is only marginal. As noted in §4.2| 
a high fraction of the clusters have significant substructure, 
but none of them shows a clear bimodal distribution in veloc- 
ity and we therefore assume that Eq. [T]is applicable to all the 
clusters in the sample. 

would be <g 1%, and still negligible over the statistical uncertainty in the 
velocity dispersion. 



The radius r2ooe is also listed for each cluster in Table [3] 
These have been calculated using M200C an d assuming spheri- 
cal clusters (i.e., M^ c = 200p c x 47rr| 0()( ,/3). 

4.2. Substructure 

It is becoming widely accepted that substructure is a com- 
mon feature of galaxy clusters, and that its presence (or lack 
thereof) is related to the degree of relaxation and hen ce the va- 
lidity of the hydrostatic equilibrium hypothesis (e.g., Battaglia 
et al.||2011| and references therein). While X-ray observa- 



tions can reveal the presence of substructure in the plane 
of the sky, velocity information can reveal substructure in 
the radial direction. From X-ray observations ove r a wide 
range in masses at z < 0.3, |Schuecker et al7| ( |200 1 ) > find that 
(52 ± 7)% of galaxy clusters present significant substructure. 
Girard fet al.| ( 1997 1 find that out of 44 optically selected local 
(z < 0.15) clusters, 15 (38%) sho w significant signs o f sub- 
structure based on their dynamics. Girardi et al. ( 1997| l argue, 
on the other hand, that substructure found in clusters that show 
a unimodal velocity distribution (i.e., where the substructure 
is not of comparable size to the cluster itself) does not influ- 
ence the velocity dispersion (hence mass) measurements. 

In general, a non-negligible fraction of the galaxy clusters 
in a sample will have biased mass measurements due to sub- 
structure. These results highlight the need for a correct es- 
timation of the degree to which galaxy clusters seem to be 
relaxed or in the process of merging. 

One very basic test for substrucure involves the distribu- 
tion of measured velocities. In fact, however, none of our 
velocity histograms shows clear evidence for a bi- or multi- 
modal distribution and the velocity dispersions Sbi are con- 
sistent with Gaussian dispersions, in all cases, within la. So, 
in the following, we employ three specific tests that take ad- 
vantage of the 3-dimensional information provided by optical 
spectroscopy to assess the dynamical state of the clusters from 
a wide perspective. Table[4]summarizes the substructure anal- 
ysis. 

4.2.1. Id: BCG Peculiar Velocity 

For a cluster that is relaxed, the peculiar velocity of the 
BCG should be clo se to zero (|Quintana & Lawrie||1982| 
|Oegerle & Hill|2001| but see|Pi mbblet et al. ( 2006) for a likely 
counter example). |Oegerle & Hill| ( |2001[ ) find that the dis- 
persion of BCG peculiar velocities is ~ 160 kms -1 fo r a me- 
dian 5bi ~ 8 00 kms" 1 . Using a sample of 452 clusters, [Coziol| 
|et al.| ( |2009| l find that BCGs have a median peculiar velocity 
0.325bi and that 41% of BCGs have velocities different from 
zero at the 2er-level, but note that this number is compara- 
ble to the fraction of clusters that show signs of substructure. 
In summary, velocities consistent with zero are not necessar- 
ily expected. Dominant (D/cD) BCGs, however, are mostly 
found in the low peculiar velocity regime. Clusters are (pro- 
visionally) considered as disturbed if their BCG has a peculiar 
velocity diffe rent from zero at the 2er-level where, following 
Coziol et al. ( 2009| l, the fractional uncertainties are given by 



A(vp ec /S B i) = — W(Av pec ) 2 + 



>BI 



Vp ec A5BI 



(2) 



where At 



pec 



5|j/A^ a i + (AvBCG) 2 is the error in the pe- 



culiar velocity, and Avbcg is twice the cross-correlation er- 
ror estimated by RVSAO, which is a conservative correction 
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Table 3 

Dynamical Properties of ACT 2008 Clusters 



ACT Descriptor 



No, 



(kms^ 1 ) 



M 2 Mc 
(1O 14 /^M ) 



ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 
ACT-CL 



JO 102- 
J0215- 
J0232- 
J0235- 
J0237- 
J0304- 
J0330- 
J0346- 
J0438- 
J0509- 
J0521- 
J0528- 
J0546- 
J0559- 
J0616- 
J0707- 



-4915 b 

-5212 

-5257 

-5121 

-4939 

-4921 

-5227 c 

-5438 

-5419 d 

-5341 c 

-5104 f 

-5259 s 

-5345 h 

-5249' 

-5227 

-5522 



55 
64 
82 
65 
71 
71 
88 
65 
76 
24 
55 
48 
31 
18 
58 



0.87008 ±0 
0.48009 ±0 
0.55595 ±0 
0.27768 ±0 
0.33438 ±0 
0.39219±0 
0.44173 ±0 
0.52973 ±0 
0.42141 ±0. 
0.46072 ±0 
0.67549 ±0 
0.76780 ±0 
1.06628 ±0 
0.60910±0 
0.68380±0 
0.29625 ± 



00010 
00012 
00009 
00006 
00009 
00008 
00009 
00007 
00011 
00006 
00032 
00010 
00020 
00026 
00044 
00006 



1321 ± 106 
1025 ±102 

884 ± 110 
1063 ± 101 
1280 ±89 
1109 ±89 
1238 ±98 
1075 ± 74 
1324 ±105 

846 ±111 
1150± 163 

928± 111 
1082± 187 
1219± 118 
1124± 165 

832 ±82 



1789 ± 140 
1736± 173 
1438 ± 177 
2007 ±190 
2339 ± 162 
1971 ± 155 
2138± 166 
1770 ±122 
23 10 ±182 
1451 ± 189 
1744 ±245 
1337 ± 159 
1319 ±226 
1916± 184 
1699 ±244 
1561 ± 156 



16.3±3.8 
9.6±2.8 
5.9 ±2.2 

11.9±3.4 

20.0 ±4.2 
12.7±3.0 
17.1±4.0 
10.7±2.2 

21.1 ±5.0 
5.5±2.1 

12.1 ±5.1 
6.1±2.2 
8.1±4.2 
14.9 ±4.3 
11.2±4.9 
5.7 ± 1.7 



a Number of s pectroscopically confirm ed members, after app lying the selection proce dure of j |3.2| 
b "El Gordo" jlVlenanteau et al.|201 lj; SPT-CL J0102-4915 jWilliamson et al.|201l} . 
c Abell 3128 (NE) (Werner et al.|2007). 

d PLCK 0262.7^10.9 jPlanck C ollaboratio n|201 la) , SPT-CL J0438-5419 ([ Williamson et al.|201 1 
c SPT-CL J0509-5341 (Staniszewski eta' '* 

f 



SCSO J052113-510418 (Menanteau et al. 2010al, SPT-CL J0521-5104 (Vanderlinde et al. 20101. 
SPT-CL J0528-5259 I Staniszewski et al. 2009 1, SCSO J052803-525945 (Menanteau et al. 2010al. 
SPT-CL J0547-5345 I Staniszewski et al. 2009 1. 



1 SPT-CL J0559-5249 (Vanderlinde et al. 20101. 

HQuintana et al.|2000] >. 

Eight clusters meet this criterion, which will be coupled 
with similarly chosen criteria in the 2d and 3d analyses be- 
fore selecting which clusters have significant evidence for 
substructure. 

4.2.2. 2d: Projected BCG-SZE Offset 

Under the hypothesis of hydrostatic equilibrium, galaxies 
closely trace the total mass distribution in the cluster and thus 
the BCG is located at the peak of the gravitational potential. 
If the cluster is virialized, the gas should also follow the mass 
distribution. Deviation from this scenario may be quantified 
by an offset between the BCG (i.e., dark matter) and the SZE 
(i.e., gas) peak. This, of course, is sensitive to offsets pro- 
jected in the sky, unlike the preceding and following tests. 



ACT h as a beam of l'A (FWHM) at 148 GHz ( |Hincks et 
|al.]2010| and the uncertainties in the determination of the po- 
sition of each cluster are of order 10"-15". We therefore list 
the projected offset in arcsec in Table E| offsets < 15" are 
within ACT's positional uncertainty andshould therefore not 
be considered physical offsets. Column 6 of Table |4] lists the 
projected offset between the BCG and the SZE peak for each 
cluster relative to th e characteristic scale of the cluster r2ooc- 
|Lin & Mohr| ( |2004l ) find that > 80% of BCGs are o ffset from 
the peak gas emission by Ar/r2oo c < 0.2. Moreover, |Skibba et] 
[aLl ( |2011| l find that ~ 40% of BCGs do not sit at the minimum 
of the potential well in clusters. 

We choose Ar/r2oo e ~ 0.20 as the threshold between (ten- 
tatively classifi ed) relaxed and dis turbed clusters, based on 
the results of |Lin & Mohr| ( |2004[ l. In this case, only three 
clusters— ACT-CL J01024gT57~ACT-CL J0509-5341, and 
ACT-CL J0528-5259— have values over the threshold. Given 
that the chance of line-of-sight (l.o.s.) substructure should be 
the same as that of substructure in the plane of the skjj^] this 

6 In fact, the latter should be approximately twice as large, given the num- 
ber of dimensions covered by the plane of the sky and the l.o.s.. 



might be too stringe nt a limit. We note that the findings of 
Skibb a~e"t al.| ( |201 l| l argue that this might not be a very reli- 
able test for substructure, but we include it for completeness. 

We note that these three clusters have offsets on the order 
of an arcminute, far beyond uncertainties in the ACT SZE 
centroids and therefore qualify as physical offsets. 

4.2.3. 3d: DSTest 

By studying a large sample of statistical tests for substruc- 
ture in galaxy clusters, Pinkne y" ;t al.|([T9 96) have shown that 
the DS test (Dressier & Shectman 1988) is the most sensi- 
tive test when used individually. The test has the ability not 
only to detect the presence of substructure, but also to locate 
the latter in projected space (in the ideal cases of substructure 
not overlapping with the main system neither in velocity nor 
in projected space) and is based in the detection of localized 
subgroups of galaxies that deviate from the global distribution 
of velocities by use of the parameter A = S,(5,, where 



5 



ocal 



[(Vi-vf + io-i-a) 2 ]' 



(3) 



is computed for each cluster member, where v, and cr, are the 
mean and standard deviation of the velocity distribution of the 
TViocai members closest to the ith member, and v and a are the 
mean and standard deviation of the velocity distribution of all 
the cluster members. The significance level (s.l.) of the test is 
obtained by shuffling the velocities of each galaxy via a boot- 
strap resampling technique with 5000 iterations. Although the 
common use is that JVi OC ai = y/N gSi u in this work A is calcu- 
lated for/Viocai ranging from 5 to 12. The uncertainties in the 
s.l. are given by the second-maximum and second-minimum 
s.l. for each cluster when varying M 0C ai (i- e -> they correspond 
to ~ 75%-level uncertainties), and the central value is given 
by the median. A large uncertainty (i.e., dependence on Mocai) 
might also be indicative of substructure, but we do not include 
this in the analysis. 



8 



C. Sifon et al. 



Table 4 

Substructure in ACT 2008 Clusters 



Cluster 




z 


|Vpec| a 


|Vpec|/S B I 


Ar b 


Ar/r 2 oo£- 


s.l. (DS) 


Disturbed? d 








(kms -1 ) 




(arcsec) 










ACT-CL JO 102- 


-4915 e 


0.870 


10± 169 


0.01±0.13 


68 


0.30 


0.48-o.n 


010 


Yes 


ACT-CL J0215- 


-5212 


0.480 


1 171 ± 153 


1.14±0.19 


33 


0.12 


0.02-™? 


101 


Yes 


ACT-CL J0232- 


-5257 


0.556 


37 ± 129 


0.04 ±0.14 


35 


0.15 


11+ 011 


000 


No 


ACT-CL J0235- 


-5121 


0.278 


138 ± 137 


0.13±0.13 


44 


0.09 


04+ 01 
uu -0.03 


001 


Yes 


ACT-CL J0237-4939 


0.334 


261 ± 174 


0.20±0.14 


78 


0.16 


<0.01 


001 


Yes 


ACT-CL J0304- 


-4921 


0.392 


151 ± 157 


0.14±0.14 


22 


0.06 


04 +0 09 


001 


No 


ACT-CL J0330- 


-5227 


0.442 


424 ± 167 


0.34±0.14 


44 


0.12 


02 l+0.27 
u - zl -0.02 


100 


No 


ACT-CL J0346- 


-5438 


0.530 


263 ± 125 


0.24 ±0.12 


16 


0.06 


23+ 05 

u, ^ J -O.07 


100 


No 


ACT-CL J0438- 


-5419 


0.421 


392 ± 172 


0.30±0.13 


10 


0.02 


0.03*g 


101 


Yes 


ACT-CL J0509- 


-5341 


0.461 


361 ± 134 


0.42 ±0.17 


114 


0.46 


0.08™ 


[10 


Yes 


ACT-CL J0521- 


-5104 f 


0.676 


440 ±292 


0.37 ±0.25 


37 


0.15 




00- 


No? 


ACT-CL J0528- 


-5259 


0.768 


144 ± 177 


0.16±0.19 


50 


0.28 


30+ 07 

u,JU -O.02 


010 


No 


ACT-CL J0546- 


-5345 


1.066 


541 ± 163 


0.50±0.17 


20 


0.13 


02+ 004 

UiUZ -O.02 


101 


Yes 


ACT-CL J0559- 


-5249 


0.609 


233 ±241 


0.19±0.20 


9 


0.03 


13+ 013 


000 


No 


ACT-CL J0616- 


-5227 f 


0.684 


685 ± 268 


0.61 ±0.25 


29 


0.12 




10- 


Yes? 


ACT-CL J0707- 


-5522 


0.296 


402 ± 140 


0.48±0.18 


19 


0.05 


34+0 04 


100 


No 



Note. — Redshifts are listed for reference. 
il Peculiar velocity of the BCG in the cluster rest-frame (see j |4.2.1) . The uncertainties consider the error on the 

BCG redshift as twice that given by RVSAO. 

b Offset between the BCG and the SZ peak as found in the Y^ooc analysis (see §^ 2. 1 14 2~2lfo r details). 

c Significance level of the DS test. Uncertainties are computed at the 75% level (see j |4.2l3] for details). 

d Each ordered number represents one of the tests listed in the table: "1" means the tests shows evidence for 

substructure and "0" means it does not. 

c This cluster is classified as "disturbed" based on the results of Menanteau et al. 1 201 1 ). See text for details. 
f There are too few members observed for the DS test to be reliable. The classification is left as a tentative one, and 



these clusters are excluded from the analysis of j |5.3| (see |4.2.4) . 



As shown by Pinkney et al. (T996JI, the false positive rate 
for the DS test is < 1%, < 4% ancT"9% for a s.l. of 1%, 5% 
and 10% respectively, for member samples as large as ours. 
The threshold for substructure detection is set therefore at 5% 
s.l. within uncertainties where, given the size of our sample, 
we expect no false detections. 

4.2.4. Substructure Results 

Clusters have been identified as merging systems if they 
meet at least two of the three conditions explained above, or 
if they have a s.l. of the DS test strictly below 5% within un- 
certainties. Although the second of the three conditions de- 
pends on the projected spatial distribution, it is clear that this 
analysis is biased towards l.o.s. substructure. 

ACT-CL J0102-4915 ("El Gordo") is a special case, as it 
does not show evidence for merging from the dynamical in- 
formation alone. However, both the spatial galaxy distribution 
and X-ray surface brightness distribution reveal that this is a 
very complex system where two massiv e clusters are in teract- 
ing close to the plane of the sky (Menantea u et al.|201 l| l. 

On the other hand, ACT-CL J06 16-5227 is tentatively con- 
sidered as a merging cluster given the high peculiar velocity 
of the BCG, but the DS test was not performed for this cluster 
given the low number of members. The latter note also applies 
to ACT-CL J0521-5104, although this cluster is tentatively 
considered relaxed. Thes e two clusters have been excluded 
from the analysis of §5.3| 

The last column of Table |4] states whether a cluster is con- 
sidered to be relaxed ("No") or disturbed ("Yes") while the 
previous column lists whether each cluster shows ("1") or 



does not show ("0") signs of substructure in each of the tests, 
as defined above. Combining the three criteria used, 7 out of 
14 clusters show signs of merger activity (or 8 of 16, if we in- 
clude ACT-CL J0616-5227 and ACT-CL J0521-5104). This 
number is consistent with previou s optical and X-ray stud- 
ies of local clusters (e.g., Girardi et al.|1997] Schuecker et al. 
2001 ) and is also consistent with the X-ray follow-up of SPT 



20( 

sz: 



E-detected clusters by And ersson et al.| ( f2Q 1 1 ) >. They find 
that 9 out of 15 SZE-selected clusters show signs of substruc- 
ture based purely on X-ray morphology. 

4.3. The Influence of Emission-Line Galaxies 

Clusters of galaxies are mostly populated by passive galax- 
ies. Late-type galaxies are preferentially found in the outskirts 
of clusters and associated with infalling groups. They there- 
fore tend to show a different velocity distribution (Biviano & 



Katgert|2004| . |Girardi et al.| ( [T996| l find that 29% (53%) out 
of a sample of 17 nearby clusters show differences in the ve 



sample of 17 nearby 
locity dispersion and 24% (47%) in the mean velocity, at the 
2a (la) level. Simulations also show that, where blue galax- 
ies are found (i.e., outside the core), the y tend to have a higher 
velocity dispersion than red galaxies (Springel et al. 2001b). 
The way blue galaxies are distributed in the cluster (both in 
space and i n velocity) depends, how ever, on the history of 
each cluster (Biviano & Katgert 2004). The issue is complex; 
for example |Aguerri & S anchez-Janssen]( |2010| l find no dif- 
ference in the fraction of blue galaxies between relaxed and 
disturbed clusters. 

Although our target selection procedure should not be 
strongly biased against emission-line galaxies, the observa- 
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Figure 3. The top two panels show, for each cluster indicated on the hori- 
zontal axis, the ratios of dynamical masses ( top) and cluster redshifts (middle) 
when only the absorption-line ("red") galaxies or all galaxies are used for the 
analysis. Error bars are given by AM a n/Mjj and Azm/Zaii, respectively. The 
dashed line in each panel marks a ratio of unity. The bottom panel shows the 
observed fraction of galaxies with emission lines ("blue"). Cluster names 
have been shortened for clarity. (See the text for details.) 

tions have not been designed to study this effect and the spec- 
troscopic samples have emission-line fractions of < 10% in 
most cases. This number does not necessarily reflect the ac- 
tual fraction in the clusters and could be taken as a lower limit 
for it. In spite of all this, we briefly study the effect that blu^j 
galaxies might have on the results. 

Fig. [3] shows, in the top and middle panels respectively, 
the variation of the mass measurement and the cluster red- 
shift when blue (i.e., emission-line) galaxies are, and are not, 
included. The null hypothesis (i.e., no bias) corresponds to 
M m i/Md\ = 1. Within uncertainties (given by AM all /M a u), 
dynamical masses do not change when including emission- 
line galaxies. Some cluster redshifts, however, change signif- 
icantly. Four clusters have different redshifts at the 2cr-level. 
The two most extreme cases are ACT-CL J02 15-52 12, for 
which the redshifts calculated including, or not, blue galax- 
ies, are different at the 6.5a-level, and ACT-CL J0102-4915, 
at the 4.5<r-level. Note that for the red-only analysis, the blue 
galaxies are removed before the selection process (i.e., N Te ^ 
does not necessarily equal iVgai — iVbiue)- 

The redshifts, velocity dispersions and corresponding 
masses in Table [3] have been calculated using all galaxies, 
since blue galaxies do not bias our mass measurements. This, 
in turn, is cons istent with the findings of Aguerri & S anchez-| 
|Janssen| ( |20!()l i. 

5. SZE-MASS SCALING RELATIONS 



Both |Vanderlinde et al.] ( [20T0] > and Seh gal et al | (2011 
have shown that, given an accurate calibration of the SZE- 
Mass scaling relation, the inclusion of the ACT or SPT clus- 
ter samples can lead to significant improvements in cosmo- 
logical parameter uncertainties, particularly w and <7g, over 
WMAP-7-only constraints. These results have recently been 

7 Although the classification is done purely based on the spectral features 
of each galaxy (with- or without emission lines), we sometimes speak of blue 
and red, instead of emission- and absorption-line galaxies respectively, for 
convenience. 



confirmed by Benson et al. ( 201 l| l using X-ray observations. 
However, without a precise SZE-Mass scaling relation, these 
cluster samples do not provide significant improvements in 
constraining cosmological parameters. Determining a precise 
and accurate scaling relation characterizing the mass of a clus- 
ter in terms of its SZE signal is the main goal of this study. 

Observations have shown that the SZ E signal and mass of 
a cluster can be related by a power-law ( Benson et aL]|2004 
|Bonamente et al.|2008| |Melin et al.|201 1 ). While most simu- 



lations seem to con firm this ( |da Silva et al.|2 004 ; Mo tl et aL] 
|2005| Nagai 2006), others suggest that certain effects (e.g., 
AGN feedba ck) can cause de viations from a single power-law 
dependence (Battaglia et al. 2011). In this work, we restrict 
ourselves to a power-law relation between dynamical mass 
(see ^4] and Table [3) and each SZE estimator measured from 
the ACT data (see |2.1.2| and Table [1} of this form: 

M 200c = 10* [y ' 5 E(zT 2 ] B h~ 7 l M Q (4a) 
M 200c = 10 4 [y E(zT 2 ] B M Q (4b) 

M 200c = 10* [Y 2m D A (z) 2 E(zT 2/3 ] B M Q (4c) 

where Da(z) is the angular diameter distance in Mpc, M200C is 
in units of A^M© and£(z)= [£I m (1+z) 3 + Qa] 1/2 . We refer to 
B as the (logarithmic) slope of the scaling relations. The self- 
similar predictions are 1, 1, and 0.6 for the yp.ss yo and Y200C 
scalings, respectively (e.g.,|Bonamente e t al. 2008 ; Marriage 
|et al.|201 la| |Sehgal et al.|2011| i. We note that these are con- 
venient forms of parameterizing the scaling relations if one 
wants to predict the mass of a cluster using SZE observations. 

Before we proceed to estimate scaling laws we consider the 
selection biases, specifically Malmquist and Eddington, that 
can affect our st udy when fitting the scaling relations (see 
Mantz et al.|2010a} for a pedagogical description). The Ed- 
dington bias results fro m the asymmetry of the steep underly- 
ing mass function (e.g., Jenkins et al. 200Tj |Tinker et al.|2Q08| > 
which introduces a shift in mass due to the statistical fluctua- 
tions of the measurement of the mass proxy in the observable- 
mass rekdtion^lTiedynamical-mass scaling relations used 
here ( |Evrard et al.|20 08) were calculated from cosmological 
simulations that inherently take into account the proper asym- 
metric mass function and therefore are not subject to Edding- 
ton bias. Our sample, however, has been selected based on 
SZE signal, which makes it subject to Malmquist bias with 
the largest biases experienced by clusters near the detection 
threshold of the survey. In the following section we describe 
our methodology to correct for this bias. 

5.1. Malmquist Bias Correction 

The Malmquist bias is produced by the intrinsic scatter in 
the observable: clusters with mean SZE fluxes at the detection 
limit whose signals scatter up will make it into the sample, 
while those that scatter down will not. We use the simulations 



of Sehgal et al. ( 2010| > to investigate this effect i n our sam- 
ple, by measuring yo m the simulated clusters ( Hasselfield et 
|al.|2012, in prep.[ ) in the exact same way as for the observed 
clusters and subsequently applying a cut (yoTcMB > 150/j.K) 
to the simulated data that reflects approximately the detec- 
tion threshold of the observed cluster sample. This procedure 
mimics the observational situation with the exception that it 
assumes a constant noise level throughout the survey. The 
ACT sample is defined in terms of a S/N limit, although the 
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Table 5 

Best-fit parameters of scaling relations 



Relation 


Bss a 


Best-fit 
A B 


0~MY 


$F-Mwoc 


1 
1 

0.6 


20.0 ± 1.9 
18.2±0.7 
17.0± 1.6 


1.16±0.27 
0.74 ±0.11 
0.56±0.11 


0.16±0.03 
0.12±0.02 
0.11 ±0.02 



a Expected logarithmic slope from self-similar evolution. 

noise level is approximately constant except near the edges of 
the map (Marriage et al. 201 lab). Within our sample, there 
is potentially only one cluster (ACT-CL J0707-5522) whose 
Malmquist correction is not accurately described by this pro- 
cedure because it sits in a high-noise region in the maps. If 
this cluster is removed from the sample, the change in the 
Malmquist-corrected scaling laws is negligible. In practice, 
clusters within a mass range from M to M+ AM are extracted 
from the simulations and the average yo value of the extracted 
subsample is determined both with and without a detection 
threshold. The ratio of yo values represents a statistical esti- 
mate of the Malmquist bias factor for clusters within this mass 
range. The procedure is repeated for different mass ranges 
that span the ACT cluster sample. At the low mass end of 
the cluster sample the bias correction factor is ~ 0.8, while 
for clusters with M200C > 9 x 1O I4 M the correction factor is 
close to unity. A continuous smooth curve is fitted to the bias 
correction factors as a function of mass and applied individu- 
ally to each cluster's individual measured y values. To distin- 
guish the bias-corrected values hereafter, we label them with a 
superscript "corr". We apply the same bias correction factors 
to each of the different SZE estimators. This is a reasonable 
approach since the latter are all based on matched filters with 
kernels of similar scales. 

5.2. Best-fit Scaling Relations 

We use the Bivariate Correlated Errors and intrinsic Scat- 
ter (BCES) bisecto r algorithm for linear regression ( |Akritas| 
|& Ber shady 1996 ), which takes into account measurement 
errors in both axes and intrinsic scatter, to find the best-fit 
slopes and normalizations of the power-law scaling relations 
given by Eqs. [4] The results are shown in Fig. H] where the 
solid lines represent the best-fit power-laws and the shaded 
regions are the la uncertainties. Table Bllists the best-fit pa- 
rameters, where the last column lists thelog-normal intrinsic 
scatter orthogonal to the best-fit line, as introduced by Pratt et 
[aL] ( |2009] ) . Uncertainties in the intrinsic scatter are computed 
following Planck Collaboration (2011b). Differ ent symbols 
identify the dynamical state of each cluster (see 5 5.3 1. 

All three SZE estimators correlate with dynamical mass 
with Pearson's r-values of 0.63, 0.75 and 0.67 for y c Q 0IT 
and T 2 c q ( ", respectively. The fractional errors on the slopes 
are similar, ranging from 15% for y^* 1 - M^ c to 25% for 
yo°^-M2ooc, while (7 my ranges from 11% to 16% depending 
on the SZE estimator. These values are consistent with those 
found in simulations (which have some dependence on the 
input cluster physics) and are of o rder 10%- 15% for large- 
aperture integrations such a s F20OC (Nagai 2006, Yan g et aT] 
|2010l Battaglia et al.||2011[ ). We find that the intrinsic scat- 
ter of the y^ a - M200C relation is low and similar to that of 
the F 2 c q,^ - M2000 scaling law. While numerical simulations 
predict a higher dependence on gas physics and projection ef- 
fects for central estimates (e.g., Motl et al.|20 05; Sha w et aT] 
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Figure 4. Scaling relations bet ween SZE estimators and dynamical mass for 
y 0!5 from |Sehgal et al.|p0TT) with ACT 2008 data (top), the central SZE 
amplitude j'o (middle), and l^ooc. the Compton y-parameter integrated out to 
100c (bottom), the latter two including three-season ACT data. All estimators 
have been scaled as indicated in the axis labels (see Eqs.[4j and data points 
have been corrected for Malmquist bias as detailed in the text. Solid blue 
lines show the best-fit power laws, with the 1 a uncertainties marked by the 
shaded regions (see Table[5j. Different symbols identify whether each cluster 
is disturbed (triangles), relaxed (circles), or not classified (squares). Previous 
estimates of the T2ooc-^20(fc scaling relation are shown in the bottom panel 
with dashed and dot-dashed lines (see text for details). 
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Table 6 

Best-fit parameters of scaling relations for selected subsamples 



Relation 



Sample 



Disturbed 
Relaxed 

Disturbed 
Relaxed 

Disturbed 
Relaxed 



20.7±5.1 
20.5 ±6.6 



1.34 ±0.76 
1.29±0.98 



0.18±0.05 
0.17±0.05 



18.0± 1.3 0.69±0.19 0.13±0.04 

18.4±2.0 0.79±0.27 0.13±0.04 

16.8±2.3 0.50±0.16 0.14±0.04 

17.3±2.4 0.65±0.17 0.09±0.03 



Note. — There are 7 merging and 7 relaxed clusters. See text for 
details. 



2008), comparison of the former with s imulations should be 
done with care since, as stated in jj |2.1.2| yo here is not neces- 
sarily the projected central pressure of the cluster as is often 
defined in simulations, but the normalization of the A10 pro- 
file in the optimally filtered maps. For this reason, comparison 
with self-similarity is also not straightforward fory™ 11 — Mhjoc. 

5.3. Scaling Relations for Relaxed and Disturbed Clusters 

Table [6] lists the best-fit scaling relations when separating 
the sample in to re laxed and disturbed clusters according to 
Table |4] (see §4. 2} . The effect of disturbed clusters, if any, 
is similar for all SZE estimators and is apparent as a slight, 
but not significant, change in slope and an increase of <tmy- 
While errors on samples of this size are very large, we find 
that Fjo^— M200C has the largest decrease in scatter when in- 
cluding only relaxed clusters, and the largest boost for dis- 
turbed clusters. Also, and contrary to expectations, the change 
in slope is largest for the ^ooc - ^200c scaling relation. 

Of the 7 merging clusters, 5 are at low redshift (z < 0.5) 
and of the 7 relaxed clusters, 4 are at high redshift (z > 0.5). 
We have explored redshift-induced biases in the scaling rela- 
tions and find that they are dependent on the fraction of merg- 
ing clusters in each sample. In this way, the high-z sample 
is consistent with the relaxed sample and the low-z sample is 
consistent with the mergi ng sample. 

As explained in 30 neither ACT-CL J0521-5104 nor 
ACT-CL J06 16-5227 have been considered in the present 
analysis. Including these clusters in either sample does not 
change the best-fit parameters and only changes the intrinsic 
scatter by <0.05. 

5.4. Previous Results 



Rine s et IT] ( |2010[ ) were the first to present a comparison 
of SZE fluxes and masses derived from dynamical informa- 
tion, but their sample selection did not allow for the estimation 
of a scaling relation. Here, we review some SZE-Mass scal- 
ing relations derived from other observations or mass proxies. 
While we note that many authors have presented scaling rela- 
tions in different forms and using a variety of mass proxies, 
here we compare to those that have done so in the same form 
as is done here (i.e., correcting by intrinsic evolution). When 
appropriate, cylindrical F200C measurements have been con- 
verted to spherical equivalents applying the correction for the 
A10 profile. 

When comparing to other results we have converted to 
values calculated within r2ooe by assuming that M scales as 
M oc Y a using the respective profiles found by each study. 
The normalization was estimated as the ratio of the NFW pro- 
file used for the mass and the GNFW profile used for l200c- 

The bottom panel of Fig. [4] shows these scaling relations; 
those where masses were estimated from X-ray observations 



( |Andersson et al.||20TTJ |Planck C ollaborat ion] |2011b| ) and 
weak-lensing measurements ( Marrone e t al.|201 1 1 are shown 
with dashed lines, and the dash-dotted line shows th e resu lts 
from hydrodynamical simulations by Battagl ia et al.| ( [20TT| ). 

All the scaling relations are in good agreement wit h~that 
presented here. The scaling derived by Mar rone et al.|p0TT) 
has a shallower, but consistent, slope than all the other scaling 
relations. The scaling relation derived in th is work is in ex- 
cellent agreement with those by Andersson et al. (201 1 1 and 



Planck Collaboration] ( |201 lb| l, both derived with X-ray mass 
estimates of SZE-selected clusters. 

6. DISCUSSION 

6.1. Individual Clusters 

In this section we list clusters with notable features, includ- 
ing comparison of dynamical masses presented here with pre- 
vious estimates, where available]^] The respective original or 
alternative names can be found in Table [3] With respect to 
notes on optical featu res of these clusters, t he reader is re- 
ferred to Figs. 4-10 of M enanteau et al.| ( [2010b| l, as appropri- 
ate. 

6.1.1. ACT-CL J0102-4915 "El Gordo" 

Located at z = 0.870, this cluster has the largest SZE signal 
(it is the rightmost data point in the top and bottom panels of 
Fig. |4} and is one of the most massive clusters of the sample 
according to its dynamics. This cluster looks elongated in the 
optical (in fact, it is double-peaked in the galaxy distribution; 
Menanteau et al. 201 l| l, but there are no clear signs o f line-of- 



sight substructure from the dynamical information. In Menan- 
teau et al. (2011), we show that if the cluster is divided into 



two subclusters in the process of merging (as suggested by the 
optical data), they have a mass ratio of order 2:1, with a total 
summed dynamical mass of M200C = (24 ±7) x 10 14 IijqMq, 
making this a huge merger between two already massive clus- 
ters^ 

Menantea u et al.| ( |2010b| l used a multi-wavelength data-set 
combining X-rays, SZE and the information provided in this 
work to estimate the cluster mass using several mass proxies 
which are in statistical agreement, with a combined mass es- 
timate of M 2 ooa = 2 1 .6 ± 3 .2 x 10 14 M Q . The statistical er- 
ror for the combined mass is likely an underestimate gi ven the 
complex nature of this massive merging cluster (see |Menan-| 
|teau et al.|201 1| for a fuller discussion of the mass measure- 
ments in "El Gordo") 

6.1.2. ACT-CL J02 15-52 12 

As shown in Table H] this cluster appears to have substruc- 
ture with a high significance as given by the DS test. More 
noteworthy, however, is the peculiar velocity of the BCG, 
y pec = 1171 ± 153kms _1 , different from zero at > 7.5o\ This 
is the only cluster in our sample in which the velocity of the 
BCG is comparable to the velocity dispersion of the cluster, 
and the cluster where the emission-line galaxies are most dif- 
ferent from the whole population. ACT-CL J02 15-52 12 has 
a second galaxy ~ 23" away (corresponding to a projected 
distance of WO/z^kpc at z = 0.480) which is only 0.27 mag 

8 We quote the original mass estimates, given as M5001- — also with respect 
to the critical density o f the Universe — and assume a typical conversion 
factor M200C ~ 1 6M500C I Duffy et al. 2008 1 when comparing with our results. 
In the particular case of lil (jordo", masses are originally given as M200a, 
the masses within a radius containing 200 times the average density of the 
universe. For this cluster, the conversion is M200C ~ 0.86A/2(X)a- 
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fainter and has a peculiar velocity of roughly -660km s" 1 , 
and at least three more galaxies within 0.55 mag of the BCG 
(which is the brightest of all by definition, but is also the one 
nearest to the optical center of the cluster), all of which have 
comparable (~ 1000km s" 1 ) peculiar velocities. On the other 
hand, this cluster does not significantly depart from any of 
the scaling relations of Fig.[4j showing the complexity of sub- 
structure analyses. It also has the highest fraction of emission- 
line galaxies. 

6.1.3. ACT-CL J0237-4939 

Similar to the previous case, this cluster has 3 bright galax- 
ies within 2' of the BCG, which are within 1 mag of the BCG. 
In particular, the second-brightest galaxy is 65" (310^^ kpc) 
away from the BCG, is 0.49 mag fainter and has a peculiar 
velocity with respect to the cluster of ~ 1850kms _1 . All this 
argues in favor of the classification of this cluster as a dis- 
turbed system. 

6.1.4. ACT-CL J0330-5227 



As mention e d in §2.1 -1| this cluster was discovered by 
Werner et aL] (|2007])beTund Abell 3128 (z = 0.06) using 
XMM-Newton X-ray observations. The SZE measurement is 
clearly associated with the background structure while the less 
massive, foreground cluster has no significant SZE emission 
(Hincks et al. 2010). The dynamical mass estim ated here i s 
significantly higher than that derived by |Werner et al.| ( [2007| >, 
of M 5 o(v = 3.4 x 10 14 /i7 ( )M Q . They do caution, however, that 
their estimate is uncertain, as it is based on an isothermal beta- 
model for the cluster. Being located only 12' away from Abell 
3128 at z = 0.06 on the sky, this cluster is a clear illustration 
of the mass selection of SZE surveys, approximately indepen- 
dent of redshift. 

6.1.5. ACT-CL J0438-5419 

ACT-CL J0438-5419 is the only new ACT cluster also 
reported by the Planck sate llite (PLCKESZ G262.7-40.9, 

has been followed-up with 
anck Co llaboration (201 Id ) esti- 
mated a mass M 5 oo c = (6.9 ±0.7) x IO^/j^M© using a Y x -M 
scaling. This value is 1.6er lower than our dynamical esti- 
mate; in fact, this cluster is one of the most massive ones in 

our sample. 

This cluster is also reported in Williams on et al.| ( |201 l| l. 
They estimate a simulation-based SZE-estimated mass 
Msooc = (8.2 ±2.5) x IO^/z^Mq, consistent with both our 
mass estimates within ler. The XMM-Newton value is sig- 
nificantly lower than the other estimates, perhaps a result of 
the disturbed state of the cluster, as we propose here. 

6.1.6. ACT-CL J0509-5341 

This was one of the first clusters discovered by SPT 
( |Staniszewski et aT7||2009[l and the first mass measurements 
were reported by |Menanteau & Hug hes (2009). Vand erlinde| 
et al. ( 2010| > reported a simulation-based SZE est imate of 
the mass of M 5 oo e = (4.3 ± 1.1) x 10 14 /i^M Q and lAnders- 



Planck ToTlaboration|[2^Tla[). It 
XMM-Newton, with which |Planc]< 



son et al. (2011 1 estimated an X-ray F^-derived mass from 



M 5 ooc = (5.4±0.6) x 10 Ii^Mq. All previous values are in 
agreement with our estimate. Cons istent with our substruc- 
ture analysis, Andersson e t al.| ( 2"0TT| found that this cluster is 
a disturbed system based solely on X-ray morphology. 



6. 1 .7. ACT-CL J0521-5104 

This cluster is not reported in Marriage et al.| ( |201 la[ ), 
because it was not a S/N>3 detection in the initial analy- 
sis. However, more recent analyses including data from 3 
years of observations show that this cluster is now detected 
at 4.5cr, and it i s therefore included in this study. Vander- 
|linde et aX] ( |2010| l report an SZE-estimated mass of M500C = 
(2.97 ±0.89) x 10 14 /z^M Q , significantly lower than the dy- 
namical mass reported here. 

6.1.8. ACT-CL J0528-5259 



This cluster was also reported by Staniszewski et al. ( 2009 ) 
and characterized optically by Menanteau & Hughes (2009 



Its SZE-estimated mass i s M 50 o e = (2.9 ± 0.9) x IO^H^Mq 
(Williamson et al. 201 1|> and its X-ray-estimated m ass is 
M 500c = (3.0 ±0.9) xW A hr 1 lM Q (jAndersson et al.||201l|). 



These values are consistent with our dynamical estimate. Also 
consistent with the present finding, Anders son et al.| pOTTJ 
found that the X-ray morphology shows this cluster to be re- 
laxed. 

6. 1 .9. ACT-CL J0546-5345 

This is the highest-redsh ift cluster of the sample, at z = 
1 .066. Brodwi n et al.| ( |2010| l first presented a spectroscopic 
study of this cluster based on 1 8 cluster members, which have 
been included in this study, including the three emission-line 
galaxies not used for their mass measurement. We included 
Brodwin et al. ( 2010| Ts galaxies in our spectroscopic catalog 



and applied the cluster membership algorithm ( §3.2| i which 
resulted in 48 members in total. Our mass estimate, calculated 
now with three times as many galaxies, is consistent both with 
their dynamical estimate and their best estimate, combining 
X-ray, SZE and dynamical information, which corresponds to 
Mzoor = (7.9±0.9)x 1O 14 /^M . 

6.1.10. ACT-CL J05 59-5249 



This cluster was al so detected by SPT and reported in |Van-| 
|derlinde et al. | ( |20T0] > . They report a simulation-based SZE- 
derived mass M 5()0c = (5.3 ± 1.2) x 1O I4 ^M , while 



An- 



dersson et al.| ( |2011[ > estimate an X-ray F^-derived mass of 
Msoqc = (6.4 ±0.5) x l0 u hj{ j M Q . Both these estimates are 
consistent with each other, and combined are consistent with 
the lower limit of our dynamical estimate. The ACT SZE sig- 
nal is consi stent with the dynamica l mass (cf. Fig. [4]). Using 
X-ray data, Anderss on et al.| ( |201 1) suggest that this cluster is 
in the process of merging but our substructure analysis finds 
no evidence for substructure. These two results are not neces- 
sarily in contradiction since X-ray morphology and dynamical 
information are sensitive to substructure with different orien- 
tations. 

6.1.11. ACT-CL J06 16-5227 



The o ptical imaging of this cluster by [Menanteau et al.| 
(2010b) was sufficient to provide confirmation but shallower 
than required to secure an adequate galaxy catalog for spec- 
troscopic targeting. Out of 73 spectra obtained, only 18 
are cluster members. Another 6 are foreground/background 
galaxies. The remaining are all late-type (mostly M) stars, 
which have similar colors to the cluster members. Both the 
SZ signal and the X-rays argue in favor of this being a mas- 
sive cluster, supporting the dynamical estimate. This cluster 
readily shows the power of the SZE to reveal high-redshift 
clusters over optical selection techniques. 
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6.2. The SCS clusters 

Of the 19 clusters observed during this program (see Ta- 
ble [2]), only the 16 listed in Table [3] were detected by ACT 
The other 3 clusters are listed in Table [7] These 3 clusters 
(hereafter "the SCS clusters") were discovered optically in the 
SCS and were included in the spectroscopic progra m because 
of their high optical richness (Menantea u et al.|201 0a), along 
with ACT-CL J0521-5104. While their redshifts are consis- 
tent with the ACT sample, all three SCS clusters are in agree- 
ment with having dynamical masses low enough to be below 
the detection threshold of ACT, and provide additional evi- 
dence that the present sample is approximately mass-limited 
in nature. 



this is in broad agree ment with numerical simulations (e.g., 
Battag lia et al.|20TT) , and in excellent agreem ent with scal- 
ing relations derived with SPT and Planck data ( Andersson et 
|aTT1|20TT] |Planck Collaboration ||2011b| l. We find that the in 
trinsic scatter of the yQ OTr -M2ooc relation is similar to that of 

In summary, the first sample of ACT spectroscopically 
followed-up SZE-selected clusters has provided results that 
agree with the expectations for the first-generation of SZE 
surveys. Scaling relations derived from this sample also agree 
with the expectations, increasing our confidence in the utility 
of the SZE for detecting and characterizing massive clusters at 
all redshifts and thus constraining cosmological parameters. 



7. CONCLUSIONS 

We have conducted a large spectroscopic follow-up pro- 
gram of clusters of galaxies discovered via the Sun yaev- 
ZeFdovich effect by ACT in its southern sky survey (Mar- 
riage et al.|201 la Menanteau et al. 2010b). We used 89 hours 
of multi-object spectroscopic observations divided between 
FORS2 at VET and GMOS at Gemini-South. With a few 
(3-4) hours of observation per cluster, we have been able to 
confirm an average 65 members per cluster, which allowed 
us to: 1) obtain robust redshifts for each cluster; 2) measure 
velocity dispersions with errors < 10%, which translates to 
uncertainties of < 30% in mass estimates; and 3) determine 
the dynamical state of the clusters. 

The cluster sample spans a redshift range 0.28 < z < 1.07, 
with a median redshift z = 0.50. Careful examination of possi- 
ble substructure shows that ~ 50% of the clusters in the ACT 
sample contain significant substructure, consiste nt with the 
X-ray study of SPT SZE-selected clusters (And ersson et al.| 



|20TT 



and with optically- and X-ray-selected local clusters 
We find that the presence of emission-line galaxies within 
clusters, which could be associated with infalling groups, 
does not significantly modify the mass estimates. For this rea- 
son, emission-line galaxies have been included as members in 
the final samples. 

Dynamical masses have b een estimated from the radial ve- 
locity dispersions using the Evrar cTet al.| (2008) simulation- 
based (T-M20DC scaling relation. These clusters have masses 
5 < M200C 22 in units of 10 14 /i™M m , with a median mass 
~llxl0 14 /z^M Q 



7O J "0' 

in agreement with the mass distribution of 



the ACT sam ple as simulated by Marriage et al. ( |201 la} and 
estimated by |Menanteau et aL] ( |2010b[ ). These clusters rank 
therefore among the most massive clusters in the Universe. 

The scaling relation between dynamical mass and SZE sig- 
nal has been studied using three estimators of the SZE: the 
single-season Compton signal within a 0f5-width pixel, jors, 
the central SZE amplitude, yo, and the Compton signal inte- 
grated within r2ooc> ^200c- In order to derive unbiased scaling 
relations a simulation-based Malmquist bias correction has 
been applied to the data, and the scaling relations include in- 
trinsic evolution with redshift. 

These scaling relations are summarized in Table [5] and rep- 
resent the main result of this work. The intrinsic scatter in 
these relati ons is con sistent wi th that pred icted by simulations 
(e.g., |Motl et al.|[2005| [Reid & Sperge1J|2006] l. We find that 
F200C ls our most robust mass proxy, with a log-normal intrin- 
sic scatter of 11 ±2%. In all three SZE estimators, disturbed 
clusters do not bias the scaling relations but might slightly 
boost the intrinsic scatter (a ~ 2-5% effect), which is lowest 
(9 ± 3%) for the Yf°™-M 2 ooc relation for relaxed clusters. All 
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APPENDIX 
ELECTRONIC DATA 

Table [8]lists the properties of the BCGs for each of the ACT clusters (see Table [3}. This table is an excerpt from the full table 
available online (from which the BCGs are shown for convenience), which contains all cluster members for the 16 ACT clusters. 
It is given for guidance in its form and content. Column 1 lists the adopted identification, based on the J2000.0 position of each 
galaxy and using the initials of the first three authors of this paper to identify the catalog. Columns 2 and 3 list the positions of 
the galaxies. Column 4 lists the magnitude in the i band; column 5 lists the cross-correlation redshifts and their associated errors 
as given by RVSAO. Column 6 lists the cross-correlation S/N r cc and column 7 lists the main spectral features of each galaxy. 
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Table 8 

Redshifts and magnitudes for the BCGs of the 16 ACT clusters 





Identification 




(hh:mm:ss) 


(ddinrniiss) 


lllj 








F cc 


Main spectral 
features 


SMH 


J010257.7- 


-491619 


.2 


01:02:57.74 


-49:16:19.2 


19.186 





.87014 ±0 


,00030 


3.39 


Ca-II K,H; [Oil] 


SMH 


rm i ?n i 
JUZijlZ.J- 


en tic 
OZlZZj 


.3 


02:15:12.26 


-52:12:25.3 


18.678 





.48587 ±0 


.00016 


3.90 


Ca-II K,H 


SMH 


J023242.8- 


-525722 


.3 


02:32:42.80 


-52:57:22.3 


18.410 





.55592 ±0 


.00014 


4.53 


Ca-II K,H 


SMH 


J023545.3- 


-512105 


.2 


02:35:45.28 


-51:21:05.2 


16.493 





.27825 ±0 


,00015 


7.18 


Ca-II K,H 


SMH 


J023701.7- 


-493810 


.0 


02:37:01.71 


-49:38:10.0 


17.582 





.33554 ±0 


.00016 


10.42 


Ca-II K,H 


SMH 


J030416.0- 


-492126 


.3 


03:04:16.04 


-49:21:26.3 


17.463 





.39289 ±0 


.00020 


9.43 


Ca-II K,H 


SMH 


J033056.8- 


-522813 


.7 


03:30:56.83 


-52:28:13.6 


17.520 





.43969 ±0 


,00019 


10.23 


Ca-II K,H 


SMH 


J034655.5- 


-543854 


.8 


03:46:55.49 


-54:38:54.8 


18.577 





.53107 ±0 


.00013 


6.16 


Ca-II K,H 


SMH 


J043817.7- 


-541920 


.7 


04:38:17.70 


-54:19:20.6 


17.470 





.41955 ±0 


.00012 


9.42 


Ca-II K,H 


SMH 


J050921.3- 


-534212 


.2 


05:09:21.38 


-53:42:12.2 


18.361 





.46257 ±0 


.00022 


7.53 


Ca-II K,H; [Oil] 


SMH 


J052114.5- 


-510418 


.5 


05:21:14.54 


-51:04:18.6 


19.060 





.67780 ±0 


.00041 


3.96 


Ca-II K,H 


SMH 


J052805.3- 


-525952 


.8 


05:28:05.30 


-52:59:52.8 


19.715 





.76695 ±0 


.00037 


6.10 


Ca-II K,H 


SMH 


J054637.6- 


-534531 


.3 


05:46:37.67 


-53:45:31.3 


21.184 


1 


.06255 ±0 


.00016 


6.47 


Ca-II K,H 


SMH 


J055943.2- 


-524927 


.1 


05:59:43.23 


-52:49:27.1 


19.103 





.61035 ±0 


.00027 


3.88 


Ca-II K,H 


SMH 


J061634.1- 


-522709 


.9 


06:16:34.05 


-52:27:09.9 


18.594 





.68765 ±0 


.00011 


6.87 


Ca-II K,H 


SMH 


J070704.7- 


-552308 


.4 


07:07:04.67 


-55:23:08.4 


16.754 





.29451 ±0 


.00019 


6.05 


Ca-II K,H 



Note. — BCGs of the 16 SZE-selected clusters presented here. This table is an extract from the full table provided 
in electronic form, which contains the data for the members of all clusters. It is provided for guidance on its form and 
content. Galaxies have been named based on their positions, and using the initials of the first three authors of this paper 
to identify the catalog. 



